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INTRODUCTION N 68- 2 3 5 85 

The well-known general space vehicle trajectory estimation problem 
consists of processing a body of tracking and other observational data 
to estimate not only the parameters of the trajectory itself, but also 
a set of other uncertain parameters which are significant in the problem. 
Such other parameters include the tracking station locations, measure- 
ment biases, and assorted astrodynamic constants. 

It is also conceivable to include in this list the velocity of 
light, since its uncertainty has a significant effect on certain types 
of measurements . Thus, in principle, it would seem reasonable to add 
the velocity of light to the list of state variables, and thereby obtain 
a better estimate not only of the trajectory parameters, but also of the 
velocity of light itself. 

The purpose of this paper is to show that the hope of improving the 
knowledge of the velocity of light in this way is an illusion, that the 
velocity of light is an unobservable parameter in trajectory estimation, 
and that its inclusion in the problem is not only inappropriate but also 
can lead to erroneous results . Although many readers of this paper may 
already know and understand this, the analysis presented here is felt to 
be useful for others who have had difficulty understanding the fact, and 
also may give some insights into other observability problems. 


OBSERVABILITY IN SEQUENTIAL ESTIMATION 


Observability generally means the ability of a measurement or obser- 
vation scheme to "see" and separate out in some sense al 1 the various 
states of a multivariable system. Thus, ordinarily, observability is 
regarded as a property of the measurement scheme itself. However, there 
is more to observability than this in the case of est ima tion based on 
noisy observations. This is particularly clear for sequential estima- 
tion, as illustrated by the well-known equations for sequential minimum- 
variance processing of data in a linear system: 
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Here, x and x* are the estimates of the state vector x, respectively, 
before and after incorporating the information contained in the scalar 
observation, y = Hx + e, where e is a random observation error with 
zero mean and variance q; P is the covariance matrix of the error in 
estimate, P = E(x - x) (x - x) T . 


Observability has to do with the information content of the obser- 
vation relative to the state x. Note that if an element of the PH 1 
vector is zero, then the estimate of the corresponding element of the 
state remains unchanged, regardless of the size of the residual, (y - Hx) • 
Also, the corresponding row and col umn of AP = (PH^ 1 ) (HP)/(HPH^ + q) will 
be zero , so that no change occurs in these elements of the covariance 
matrix. For the purposes of this paper, a state element is said to be 
unobservable if the observations do not change the estimate of that ele- 
ment. Thus, so far as the observation y = Hx + e is concerned, the 
ith element of the vector PH^ is zero . 


Observability thus depends not only upon the H row vector which 
characterizes the observation, but also upon P. It is, therefore, 
crucial that the correct P matrix be used in considerin g questions of 
observability in sequential estimation. In particular, one must be con- 
cerned with the initial P matrix for a given problem, and be sure that 

it properly summarizes the a priori knowledge of the state x as it is 
supposed to. This, it will be noted, is in contrast to statements which 

sometimes have been made to the effect that the initial P matrix is 

not too material, or can be specified somewhat arbitrarily. 


SPECIFICATIONS OF A PRIORI COVARIANCE MATRIX 


For a typical trajectory estimation problem, the state elements are 
vehicle position and velocity (six elements), tracking station location 
(three elements per station), earth radius, earth-moon distance, earth 
mass, moon mass, velocity of light, etc. The initial covariance matrix 
summarizing the a priori distribution of error in the est ima tes of all 
these variables can be constructed from independent measurements of each 
of the variables . A typical example is the set of observations used 
initially to determine the tracking station locations . Assuming that 
surveying techniques were used, which employ (among other things) dis- 
tance measuring equipment (DME) , it can be argued that the uncertainty 
in the velocity of light enters in the following way. Assume that DME 
uses a two-way phase measurement. 
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<P = (3) 

where <p is the phase shift of an electromagnet ic wave of frequency f 
traversing the distance 2R, and c is the velocity of propagation* 

If the state variables are taken to be R and c, then the linearized 
form of equation ( 3 ) is 


Acp = — (4) 

c c c ' 

where AR and Ac/c may be defined as the state variables of the linear- 
ized system. In terms of the measurement of AR, equation (4) may be 
written as 


AR M = r~ ^ - R — 

w 4jtf c 


(5) 


This linearized observation equation implies an H matrix (with respect 
to the variables AR, Ac/c), which is 


H = [1 -R] 


( 6 ) 


The analysis of all prior observations on the state space is, of 
course, more complicated than the simple case described above, inasmuch 
as all observations must be expressed in terms of a common reference 
system. However, it can be shown that as long as the totality of the 
prior observations "span" the state space (that is, P Q is nonsingular), 
they can be summarized by a set of equivalent measurements having H 
matrices of the same form as ( 6 ), namely, they are direct with respect 
to all the state variables except Ac/c, and indirect with respect to 
the latter. For the three components of tracking station location, the 
complete H matrix (considering, for the moment, only those variables 
which directly affect the surveying measurements) is 


Hsta 


1 

0 


0 0 ~ x sta 

1 0 -7sta 


(7) 


001 -Zgta 


For the various astrodynamic constants considered a s imil ar analy- 
sis can be applied, the details of which are not shown here, which 
results in similar equivalent H matrices. For the radius of the earth 
and the earth-moon distance, for instance, we obtain 


Hr E = [1 

h Rem = t 1 
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and for the masses of the earth and of the moon, 

H PE * U -3he1 
~ ^1 ~3Mm^ 


( 9 ) 


(in equations (8) and ( 9 ) zero elements in H have been omitted for 
simplicity.) Equations ( 9 ) indicate that the sensitivity of the equiva- 
lent measurements of the masses to the velocity of light uncertainty 
are minus three times the respective p's. The H matrix for the 
direct a priori velocity of light measurements is simply 


H c = [0 0 ... 0 lj 


( 10 ) 


since it is not dependent on any of the other state variables. 

In regard to the a priori observations of the spacecraft trajec- 
tory parameters, a simple analysis is not possible, so we will use here 
a heuristic argument. We will assume that the trajectory parameters 
used are the components of the instantaneous position and velocity vec- 
tors, and that the time of injection is t D . The a priori (i.e., before 
t Q ) observations include all measurements made at the launch site, and 
tracking and other measurements during the launch. The totality of 
these observations can be summarized in terms of equivalent direct mea- 
surements of x v , y v , z v , x v , y v , and z v , with sensitivity to the 
velocity of light uncertainty Ac/c, just as in the case with the pre- 
viously considered variables. Analysis to determine the Ac/c sensi- 
tivity coefficients is complicated, but it is reasonable to expect that 
they would turn out to be precisely the coefficients which would result 
if real direct measurements with DME could be, and were, made. As has 
been shown, these coefficients are the negatives of the respective state 
variables. Thus, we assume that the equivalent H matrix for this part 
of the state is 



again excluding all zero columns. 

Collecting all the equivalent measurement matrices, equations ( 7 ) 
to (ll), we obtain a total a priori H matrix of the form 



( 12 ) 
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■where the last column is the set of sensitivity coefficients for the 
Ac/ c uncertainty. Associated with the set of equivalent measurements 
there is an equivalent measurement error covariance matrix, Qq, which, 
partitioned in the same manner as H Q , can he represented as 


Qo = 


Qi ] o 
P ! !c_ 


(13) 


Here, the zeroes mean that it is assumed that the random errors in the 
direct measurement of c are uncorrelated with the random errors in 
the other measurements. 


With the a priori H 0 and Qo matrices, it is a simple matter to 
form the a priori estimation covariance matrix, if a least-squares 
reduction of the a priori data is assumed: 


HqVo 1 
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VELOCITY OF LIGHT OBSERVABILITY IN TRAJECTORY ESTIMATION 


Having a representation of the a priori P Q matrix, we now consider 
the processing of a scalar observation characterized by an H matrix 

H = [H 2 j b] (15) 


where b is the (scalar) sensitivity of the observation to the 
uncertainty. It is seen that the P 0 H^ vector is given by 

m (Qi + qcaa T )H 2 T - qc&l 

P 0 H = 

-q c (a T H 2 - b) 

where the lower element is the one associated with the Ac/c variable. 
As previously stated, Ac/c is unobservable if and only if the last 
element of P q hT is zero, which requires that 

b = H 2 a (17) 


Ac/c 

(16) 


The question in considering the observability of Ac/c in trajec- 
tory est ima tion is whether or not the observations made of the vehicle 
have the property (17) • But before taking up this matter, it is 
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necessary to consider the fact that in sequential trajectory estimation 
where the epoch is current time, the dynamics of the problem camse P Q 
to change. In P Q , as represented by equation (l4) , it is noted that 
the last column (except one element) is proportional to the a vector. 
Now, since the linearized system (state) equation is of the form 

x = Fx + Gu (1 8 ) 

the P matrix, in the absence of observations, varies with time accord- 
ing to the differential equation: 


P = FP + PF T + GRG T (19) 

T 

where R = E[uu ] . From equation ( 19 ) it is easy to show that in the 
present problem the last column of P, which we may call p c (i.e., 

P = [Pi j p c ]) , obeys the differential equation 

P c = FP C (20) 


which is the same as the system equation (l8) except that there is no 
random forcing function equivalent to u. Thus, where p c is initi all y 


P c (-t 0 ) 



( 21 ) 


at a later time it would become 


p c (t) = f(t,t 0 )p c (t 0 ) 


( 22 ) 


where <b is the transition matrix of the system represented by equa- 
tion (l8). In partitioned form, 

®(t>t 0 ) = 

and 


P c (t) = q c 


Now, the element of the PH^ vector which affects the estimate of c 
is the scalar product of the p c and vectors: 

( PhT ^/c = H Pc 

= qct-H 2 $xa + b] 


-4>ja 


(24) 


$1 1 0 


(23) 


( 25 ) 
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Considering that $1 is the matrix of partial derivations of the 6 tate 
vector at t with respect to the state vector at t 0 > 


®i 


8x(t) 

^x(t Q ) 


(26) 


and H2 and b are, respectively, the partials of the measurement quan- 
tity M with respect to the state variables and Ac/c, at time t, 

■> 

,T 

I dM l+.\ 

h 2 


'dM(t) ] 

! 


b = 


equation (25) can be written as 
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which, "by factoring out the scalar dM(t)/dM(t 0 ) , becomes 
(™ T )/sc/c * 1c 
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(29) 


It is seen* therefore , that if the unohservahility condition - namely, 
the quantity in brackets in (29) is zero - at any time (after all, t Q 
is arbitrary), then it will also be satisfied at any later time. That 
is to say, the unobservability condition is time invariant, and only 
the situation at t Q need be considered to establish whether or not c 
is observable . 


Now let us consider the properties of certain types of observations. 
First, for range measurements from an earth-based tracking station, where 
range is given by 

R = [ ( x v ~ x sta) + (yv “ ^sta) + ( z v “ z sta) 1 (30 ) 


it is easy to show that the nonzero elements of the H matrix (the 
partials of range measurement with respect to all relevant state vari- 
ables) are as given in Table I. 


TABLE I 


Variable 

X V 

yv 

z v 

x sta 

ysta 

z sta 

Ac/c 

Pc/ 1 ^ 

x v 

Vv 

Zy 

x sta 

ysta 

z sta 

1 

Hr 

Xv - Xsta 

yv - ysta 

zv - Zsta 

Xv - xsta 

yv - ysta 

zy - zsta 

-R 

R 

R 

R 

R 

R 

R 
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In the table the elements of the p c vector, normalized by q c , are 
also given. The last element of Hr is -R, which is b for the 
range measurement . It is evident here that the unobservability condi- 
tion is satisfied, i.e., H 2 a = b. Note that this is the same as saying 
that the dot product of the Hr and p c vectors is zero. 

Next, consider the case of range-rate measurements. Range rate is 

R = [(xv ~ ^sta) + (yv ~ ysta) + ( z v ~ z sta) 3 ^ (31 ) 

The partials which make up Hr are more complex than those of Hr. 
Although the variables x s ta> ysta > z s ta need not be introduced as 
part of the state due to their functional dependence on x s -£ a ; ysta; 
z sta> i-t 16 convenient to do so here. The partials are given in Table II, 
together with the corresponding elements of p c / q c . 


TABLE II 


Variables 



rT 

*sta 

-T 

Rsta 

Ac/ c 

Pc/^c 



r t 

K sta 

-T 

R sta 

1 

h r 

/■*■ ' ~\^ 
/R HR\ 

\R R V 

(D T 

/- -\ T 
JR. RR ) 

\R " RV 

4 ) 

-R 


In Table II, vector notation has been used for brevity. With a little 
algebra it can be seen that HrP^ = 0 , although this is not so obvious 
as in the case of Hr. Thus, Ac/c is unobservable with R measurements. 

The results indicated by the analysis given here have been verified 
by means of a computer simulation of range/range-rate tracking of a space 
vehicle, which shows that range/range-rate tracking does not give an 
improvement in the knowledge of c provided that the correct initial P 
matrix is used. Results have also been obtained which show that if the 
initial correlations between Ac/c and the other state variables are 
selected to be zero, then an order of magnitude (or better) improvement 
is indicated in the knowledge of c in certain trajectory estimation 
situations. This indicates clearly that the specification of proper 
initial correlations is crucial. 


Now consider some angle-type measurements . One case is the measure- 
ment of the antenna pointing angle for an earth-based tracking station. 
This is equivalent in a planar problem to the measurement of the angle 
between the range vector R and the earth- vehicle vector R v : 


6 = 


cos -1 



( 32 ) 
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For a simple planar problem, the corresponding H matrix is given in 
Table III . 


TABLE III 


Variable 


y v 

x sta 

ysta 


Pc/lc 

*v 

JV 

x sta 

^sta 

B 

h 6 

(yv - ysta) 2 

(x v - x 6 ta)(yv - ysta) 

(yv - ysta) a 

(x v - XetaKyv - ysta) 

0 

R 3 sin 6 

R 3 sin 8 

R 3 ein 8 

R3 ein 8 


As is seen , the measurement of S is not affected by the Ac/ c uncer- 
tainty, that is, the b element of H 5 is zero. Nevertheless, we see 
that HgPc = 0 , and 6 observations, therefore, give no information on c. 

On-board-type angle measurements can also be considered. Two types 
are possible: (l) the direction of the line of sight from vehicle to 

earth (or other planet) , as obtained from sextant or theodolite observa- 
tions, and (2) the subtended angle of the earth (or planet). Considering 
the following three measurements , 

a = declination 

(3 = right ascension ) of earth from vehicle 

7 = half-subtense 

j 

the corresponding H matrix is given in Table IV. 


TABLE IV 


Variable 

x v 

Yv 

z v 

Re 

Ac/c 

Pc/'lc 

*v 

yv 

z v 

Re 

■ 


XyZy 

Yv z v 

(*V + yv ) 1/2 

0 

0 

R?(x? + y?) 1/a 

R^(^ + 


Hp 

Yv 

4 + yf 

*v 

(xv + yv) 

0 

0 

0 

h 7 

Re*v 

Rv 

R E Yv 

Rv 

R E Z V 

R 2 

1 

0 


Again, it is seen that such observations contain no information relative 
to the velocity of light. 
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Although the observations described here are the most common types 
considered in trajectory estimation, this is by no means an exhaustive 
list. Furthermore, no easy generalization is apparent. Hence, it is 
not possible at this point to state unequivocally that there are no 
observation types in trajectory estimation which would improve the knowl- 
edge of the velocity of light . The pattern is convincing nonetheless . 
What appears to be the case is that only direct measurements of the 
velocity of light, that is, repetitions of the laboratory experiments, 
can be expected to improve the estimate of c. 


ELIMINATION OF Ac/c FROM THE ESTIMATION PROBLEM 


From the foregoing it appears that there is no point in complicat- 
ing the trajectory estimation problem by including the velocity of light 
as a parameter to be estimated. One way of eliminating this parameter 
follows directly from the preceding analysis . It is clear that repeated 
observations of the trajectory improve the knowledge only of a subspace 
of the state space, and are equivalent to reducing the Qi observation 
error covariance matrix in equation (l4) . The limiting P matrix thus 
can be obtained by letting Qi approach zero in equation (l4): 


■min 


<lc 




( 33 ) 


Obviously, observations having H matrices of the form H = [H 2 J Hga] 
have no effect on P m-in , that is , Pmin^ = 0 • Thus , Ac/ c can be elimi- 
nated from the estimation problem by subtracting P m-in from the initial 
P 0 matrix. This leaves only zeroes in the last row and column of the 
reduced P matrix, and since these remain zero throughout the problem 
they need not be carried. At any time the P m in matrix can be added 
onto P, of course, if a representation is desired of the total 
uncertainty in the estimate of the state . 


Another way of eliminating the velocity of light from consideration 
is by using a different definition of the problem at the outset. It 
should be recognized that all of the really accurate distance measure- 
ments utilize electromagnetic radiation, and, in effect, measure not 
distance as such but rather the time it takes light to traverse the dis- 
tance . Thus , the scale of things in the universe must always be uncer- 
tain to the same degree as is the velocity of light. It should be noted, 
however, that this uncertainty really is a matter of the definition of 
the unit of length- Perhaps we should have noted in the beginning that 
the modern definition of the meter is so many wavelengths of the orange- 
red line of Krypton 86, not the platinum bar in a vault in France. In 
terms of the modern standard, then, no observations, a priori or other- 
wise, are affected by the uncertainty in the velocity of light, and the 
latter simply does not enter the picture . What this means is that in 
terms of the modern standard meter the scale of things is known more 
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accurately than in terms of the old meter, 
tainty in the calibration of the two meters 
which is the same as the uncertainty in the 


The difference is the uneer- 
relative to each other, 
velocity of light . 


CONCLUSIONS 


It appears that the velocity of light is not observable in trajec- 
tory estimation problems - or for that matter, in any other problem in 
which the velocity of light enters only indirectly. It should, there- 
fore , be eliminated from consideration in such problems . 

It has also been shown that when the velocity of light uncertainty 
is considered it is crucial that the correct a priori correlations be 
used in the estimation covariance matrix. Otherwise the results would 
show an erroneous "improvement" in the knowledge of the velocity of 
light . This fact suggests the probability that proper a priori correla- 
tions are important for other variables as well in order not to obtain 
erroneous, overly optimistic results in estimation problems. 



